First order Linear Systems 



1. Models: Two Examples 

Example 1. Farmer Jones and farmer McGregor each have a field full of 
rabbits; farmer Jones' field contains x(t) rabbits, and farmer McGregor's 
y(t) ones; t, time, is measured in months. These rabbits breed fast: they 
show a net growth rate of 0.5 rabbits per rabbit per month. The rabbits can 
also hop over the hedge between the fields; the grass is greener in Farmer 
McGregor's field, so Jones' rabbits jump over a the rate of 0.2 per month, 
whereas McGregor's jump at the rate of 0.1 per month. 

We can illustrate this with the following flow diagram: 
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Putting everything together, the equations governing x and y are: 

x = .5x — .2x + .ly = 3x + .ly 
y = .5y — .ly + 2x = .2x + Ay. 

This is an example of a first order 2x2 linear system. 
Example 2. Consider the equations: 

R' = \R + J 

r = -g* + !/ 

This is another first order 2x2 linear system. It is the result of an anal- 
ysis by the MIT Humanities Department of the plot a famous Shakespeare 
play: R denotes Romeo's love for Juliet, and / Juliet's love for Romeo. What 
does this model mean? Let us try to work backwards. 

The change in Romeo's feelings towards Juliet is mostly determined by 
how she feels about him: / is the most important term in the expression for 
R'. His own feelings have a small reinforcing effect corresponding to the 
factor of 1/4 in front of the R. 
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On the other hand, Juliet is more complex. She has a healthy self- 
awareness. If she loves Romeo, that very fact causes her to love him more: 
this is where the (3/4)/ term comes from. On the other hand, if he seems 
to love her, she gets frightened and starts to love him less aAS, hence the 
-(17/16)Rterm. 

We shall revisit this example in sections text companion matrix and ap- 
plet companion matrix, and analyze mathematically the solutions to this 
tragic model. 

2. Solving a Linear System by Elimination 

Definition. A two-by-two first order linear system of ODE's with constant 
coefficients is a collection of equations: 

x = ax + by 
y = cx + Ay, 

where a, b, c, and d are constants. 

We'll shorten this to first order linear system or even linear system. 

Remark. Say you have dependent variables X\, x„; one can define 
n x n first order linear systems (with constant coefficients) in exactly the 
same way: 

k\ = flnxi + . . . + a\ n x„ 

x n = a n \X\ + . . . + a nn x n , 

where the are constants. 

In the next few sessions, we shall develop many tools for understand- 
ing 2x2 linear systems, both analytically and qualitatively. All these tech- 
niques generalize, in a fairly straightforward fashion, to n x n systems. 
However, this goes beyond the scope of this course. The interested reader 
could for instance consult the textbook by Edwards and Penney. 

The naive way to solve a linear system of ODE's with constant coef- 
ficients is by eliminating variables, so as to change it into a single higher 
order equation, in one dependent variable. One then solves this equation 
using the techniques for constant-coefficient ODE's learned in unit 2. This 
is best illustrated with a worked example. 
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Example. Let us consider the system of equations from example 1. 

x = 0.3x + 0.1y (1) 
y = 0.2x + OAy. (2) 

Step 1. Transform the equations to get a second order ODE for x. 
Use (1) to express y in terms of x: 

y = 10x- 3x (3) 

Plug this into (2). (From equation (3) we get y = 10'x — 3x). This gives 

lOx - 3x = 2x + 4x - 1.2x lOx - 7x + x = 0. (4) 

Step 2. Solve the ODE for x. 

The characteristic equation for (4) is: 

10s 2 -7s + 1 = & s 2 -. 7s + .1=0, 

which has roots T\ = .5 and = .2. Thus we get two basic solutions, 

x 1 = e 3t and x 2 = e 2t . 

Step3. Solution for y. 

Each basic solution for x gives a corresponding solution for y, using equa- 
tion (3) 

yi = 2e- 5t and y 2 = -e 2t . 

Step 4. Using superposition we get the general solution 

x(t) = Cl e°- 5t + c 2 e 02t 
y(t) = 2 Cl e°- 5t - c 2 e 02t 

Remarks. 1. It is important to understand that the constants C\ and c 2 are 
the same for x and y; this follows from equation (3). 

2. For certain c„ there will be negatively-valued solutions; these are clearly 
not biologically significant: the model only holds for x,y > 0. 

3. We chose to eliminate y to have a second order equation in terms of x; 
we could just as well have chosen to eliminate x to get an equation in y. It 
might sometimes be computationally easier to go one way than the other; 
look out for this. 
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4. We started by solving systems by elimination because it reduces to our 
previous methods. This will not be our preferred technique. In fact, in both 
theoretical and especially numerical work it is usually preferable to go the 
opposite way and convert a higher order ODE into a system of first order 
equations and then use matrix methods. 
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